Millikelvin de Haas— van Alphen and Magnetotransport studies of Graphite 
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Recent studies of the electronic properties of graphite have produced conflicting results regard- 
ing the positions of the different carrier types within the Brillouin zone, and the possible presence 
of Dirac fermions. In this paper we report a comprehensive study of the de Haas-van Alphen, 
Shubnikov-de Haas and Hall effects in a sample of highly orientated pyrolytic graphite, at tem- 
peratures in the range 30 mK to 4 K and magnetic fields up to 12 T. The transport measurements 
confirm the Brillouin-zone locations of the different carrier types assigned by Schroeder et al.—: 
electrons are at the if-point, and holes are near the H-points. We extract the cyclotron mass and 
scattering time for both carrier types from the temperature- and magnetic-field-dependences of the 
magneto-oscillations. Our results indicate that the holes experience stronger scattering and hence 
have a lower mobility than the electrons. We utilise phase-frequency analysis and intercept analysis 
of the 1/B positions of magneto-oscillation extrema to identify the nature of the carriers in graphite, 
whether they are Dirac or normal (Schrodinger) fermions. These analyses indicate normal holes and 
electrons of indeterminate nature. 

PACS numbers: 71.20.b, 71.18.+V, 81.05.U- 



I. INTRODUCTION 

Graphite is a fascinating material whose novel elec- 
tronic properties have been extensively studied. 2 It con- 
sists of weakly-bonded layers of graphene, resulting in 
a highly anisotropic Fermi surface and semi-metallic 
properties. Recently there has been a resurgence of 
interest in graphite due to the possible occurrence of 
quasi-relativistic graphene physics in this 3D bulk mate- 
rial. There is however much discussion over whether the 
Dirac fermions found in graphene are actually present in 
graphite samples. 3 ~ 9 

Tight-binding calculations using a 2D model of 
graphite were originally utilised by Slonczewski and 
Weiss-iS (SW) in 1958. They introduced the first 3D 
model of graphite, a k • p perturbation calculation using 
the tight-binding wavefunctions of the 2D model as basis 
functions. The result was a band structure dependent 
on seven parameters (70 to 75 and A) to be determined 
experimentally. 

Table U summarises the results of subsequent studies 
of the band structure of graphite, which utilise magneto- 
oscillatory effects such as the de Haas-van Alphen 
(dHvA) and Shubnikov-de Haas (SdH) effects. Mc- 
Clure— used an interpretation of previous measurements 
of the dHvA effect in single-crystal natural graphit o 12 i 13 
to determine, or place limits on, the SW parameters. He 
derived a Fermi surface of the form of three approxi- 
mately ellipsoidal surfaces along the HKH edge of the 
Brillouin zone. This analysis was unable to specify the 
parameter 72, determining the carrier types (electrons 
or holes) at the three extremal orbits, but favoured the 
assignment of holes to the extremal orbit at K and elec- 
trons to the extremal orbits close to H . Figure Q] shows 
the Brillouin zone of graphite with the ellipsoidal Fermi 
surfaces along the HKH edges. Note that in the figure 
the carrier labelling reflects that of the current accepted 




FIG. 1. Brillouin zone of graphite. Electron and hole pock- 
ets formed by the n bands are centred along the zone edges 
HKH. 



associations from the work of Schroe der et al. —, which 
contradict the original assignments of lMcClurel 

Subsequently, Soule et al.— used angle-dependent 
SdH measurements on single-crystal natural graphite to 
demonstrate that the Fermi surfaces for both holes and 
electrons are closed. They found that the extremal cy- 
clotron orbit at the K point is larger than the orbit near 
the H point. 

Williamson et al.— were the first authors to consider 
the phase of the magneto-oscillations. They followed the 
formula of Lifshitz and Kosevich— from which the fun- 
damental harmonic of the oscillatory magnetic suscepti- 
bility can be written as 

'B 
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(1) 



where Bq is the fundamental frequency of the magneto- 
oscillations for a given carrier type, and B is the ap- 
plied magnetic field. The phase factor 7 comes from the 
Onsager-Lifshitz quantization condition and is ^ for non- 
relativistic free electrons. The offset 5 is related to the 
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TABLE I. Comparison of the conclusions of studies of the carriers' locations within the Brillouin zone and their natures in 
graphite, /fft is the frequency of the 1/B oscillations for each carrier type. The phases S and 7 are defined in the text. In the 
case of Williamson, the authors did not discuss carrier natures, but their phases are consistent with the interpretations stated. 
In some other cases authors' conclusions are not consistent with the phases measured. These are marked (*) and discussed in 
the text. 



curvature of the Fermi surface in the /c z -direction, and 
is — 4 or +i for max imum or minimum extremal cross- 
sections, respectively. IWilliamson et all showed that, ac- 
cording to the SW model, 7 takes on its usual value of | 
at the K point of the Brillouin zone, but is zero at the H 
points. However, at the positions of both the extremal 
orbits (K, and -70% of the way from K to H)±±. the 
SW model predicts 7 = \ - a resu lt with which thei r ex- 
perimental dHvA results agreed. IWilliamson et al\ also 
pointed out a fundamental difficulty in measuring the 
phase: one needs to approach the quantum limit of low 
Landau-level filling factors in order to obtain accurate re- 
sults, but in this limit the oscillations cease to be periodic 
in 1/B. 

In 1968, the magnetoreflection studies of Schrocdcr 
et al.— led to a reassessment of the SW band parame- 
ters, with the conclusion that the majority carrier types 
assigned previously were incorrect. It was shown that 
electrons occupy the orbits at the K point while holes 



occupy those near the H point (as shown in Fig. [TJ. 
This was confirmed by the magnetoresistance, Hall ef- 
fe ct, thermo power and thermal resistivity measurements 
of lWoollaml ^ 

Since the fabrication of individual graphene layers^ 
the experimental interest in graphene has increased 
tremendously. Experiments have shown that the charge 
carriers in graphene are massless, quasi-relativistic Dirac 
fermions (DFs) with a linear dispersion relation result- 
ing in an anomalous quantum Hall effect with plateaus 
at half-integer filling factorsJ^^ Recently Luk'yanchuk 
and Kopelevich— & have reported 2D-like electronic prop- 
erties in highly oriented pyrolitic graphite (HOPG) and 
have claimed also to have observed the presence of DFs. 
In 2004 Luk'yanchuk and Kopelevich— presented dHvA 
and SdH experiments on a sample of HOPG at a tem- 
perature of 2 K using magnetic fields up to 9 T. They de- 
termined the nature of the carriers, either massless DFs 
or massive Schrodinger fermions (SFs, also described in 
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TABLE II. Phases 7 and S from Eq. JTJ for different carrier 
natures: Schrodinger or Dirac fermions (SF or DF); two- or 
three-dimensional (2D or 3D). 



the literature as "normal" fermions), by two-dimensional 
phase-frequency analysis of the complex Fourier trans- 
forms of the quantum oscillations. Table |TT] details the 
values of the phases 7 and 5 from Eq. [T] expected for 
carriers of different nature (DF or SF) and dimensional- 
ity (2D or 3D). They concluded that the high-frequency 
carriers (which they assigned to be holes) are 2D DFs 
and the low-frequency carriers (which they assigned to 
be electrons) are 3D SFs. In this terminology the con- 
clusion of the previous studies of Williamson et al.— 
were that both carrier types were 3D SF, raising a con- 
troversy over the nature of t he high-frequency carriers. 
ILuk'vanchuk and Kopelevichl assign ed the carrier types 
in contradiction to ISchroeder et all and IWoollaml . 

In 2006, Luk'yanchuk and Kopelevich— detailed 
further SdH experiments under similar conditions. 
The carriers w ere a gain assig ned in contradiction to 
ISchroeder et all and IWoollaml . Plotting the extrema 
of the oscillations against Landau-level index, the total 
phase (6 — 7) for each carrier type was found using the 
equation 



Aa xx (B) ~ A (B) cos 2rr 



Bo 



(2) 



where Aa xx is the oscillatory part of the longitudinal 
conductivity and A (B) is the non-oscillatory amplitude. 
Hence the carrier nature and dimensionality was deduced 
from Table [H] The conclusion was in agreement with 
the authors' previous letter^ they found that the high- 
frequency carriers (which were assigned to be holes) are 
2D DFs and the low-frequency carriers (which were as- 
signed to be electrons) are 3D SFs. This intercept anal- 
ysis method can also be applied to dHvA oscillations 
through Eq. (P). 

Mikitik and Sharlai— identified an incorrect assertion 
by ILuk'vanchuk and Kopelevichl that there is an inher- 
ent phase shift of it in the magnetic susceptibility for 
holes compared to that of electrons. They re-analysed 
the data, concluding that both carriers are 3D SFs. 
However, the phase for the high-frequency carriers (elec- 
trons), 0.5 [2tt] would be more consistent with 2D SFs. 
They also related the value of 7 to the number of band- 
contact lines encircled by the carrier orbits in fc-space: 
7 = 2 or when an even or odd number of band-contact 
lines are encircled, respectively. Both extremal orbits 
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FIG. 2. Graphical representation of the total phases cp = 
(5 — 7) [2tt] determined by the authors listed in Table [I] (a) 
for high-frequency carriers, and (b) for low-frequency carriers. 
The phase starts at zero at the top of each circle and increases 
in the c lockwise direction. In the case of the measurements of 
Schneider et all . (N) refers to their results on natural graphite 
and (H) refers to their results on HOPG. 



in graphite enclose four band-contact lines and so one 
should expect 7 = 5 (the carriers should be SFs). 

Orlita et al. * performed optical magnetotransmission 
experiments on HOPG. These measurements are sensi- 
tive to the carriers at the H and K points, not at the 
extremal cross-sections probed by dHvA and SdH. They 
confirmed that the carriers at the H point (assumed to be 
holes) had a \JB Landau-level spectrum and are there- 
fore DFs, with some indication of a 3D nature. A linear 
Landau-level spectrum was also observed and assigned to 
electrons at K. 

In 2009, Schneider et al. 7 performed SdH experiments 
on both natural graphite and HOPG at millikelvin tem- 
peratures. They utilised the phase-frequency analysis 
technique to determine carrier nature. In HOPG their 
results were consistent with the conclusion that both car- 



riers are 2D SFs. In natural graphite, they concluded 
that both holes and electrons are 3D SFs. However, 
the total phase, (6 — 7) for the electron oscillations was 
-0.28 [2tt]; 5 = -f , and so 7 ~ -0.16 [2tt] which is 
not consistent with either carrier nature, though it is 
closer to DF. Furthermore, the total phase found for the 
holes was —0.43 [2tt]; with 6 — ~\ this gives 7 ~ —0.31 
which is also no t consistent with either carrier nature. 
ISchneider et all once again raised the issue of whether an 
intercept analysis of the type used by Luk'yanchuk and 
Kopelevich 4 can yield an accurate value of the phase. 
This remains a hotly contested issue^^ 

Figure [2] is a graphical representation of the total 
phases tp reported for the carriers in graphite by the 
authors listed in Table HJ The phases for 2D DFs and 
SFs are shown at the tops and bottoms of the phase 
circles, respectively. The phases for 3D carriers whose 
orbits are at maximal Fermi-surface cross-sections are 
also marked. For the low- frequency carriers (Fig. I2bl) . 
there is reasonable consensus that the carriers are SFs, 
with some disagreement over their dimensionality, prob- 
ably because different types of graphite (natural, py- 
rolytic and HOPG) were used in different studies, and 
these had different inter-layer coupling. However, for the 
high-frequency carriers (Fig. I2"a|) . there is no consensus 
with several authors, using different techniques, report- 
ing phases inconsistent with either carrier nature. It is 
clear that there remains a controversy over the nature 
(DF or SF) of carriers in graphite, and also some con- 
cerns over the methods used to determine these from the 
phase of dHvA and SdH oscillations. To resolve this, it 
is necessary to apply all the experimental techniques and 
analyses to the same sample. In this paper we present 
the results of dHvA, SdH and Hall experiments carried 
out on the same sample of HOPG, covering a broader 
range of magnetic fields and temperatures than previous 
studies. 



II. EXPERIMENTAL DETAILS 

The sample used for the experiments detailed in this 
paper was a piece of ZYB grade HOPG measuring 10 x 
5.5 x 0.5 mm. Magnetization measurements were made 
using a torsion-balance magnetometer shown schemati- 
cally in the inset to Fig. [3] and described in detail by 
Matthews et al.~ The graphite sample had the normal 
to its atomic planes tilted, in two separate experiments, 
at 2° and at 20° to the applied magnetic field. The 
magnetometer rotor was balanced using a piece of semi- 
insulating GaAs of similar mass. Subsequently, the Hall 
and SdH measurements were performed on the same sam- 
ple, contacted using silver paint, in six-terminal Hall-bar 
geometry. Standard low-frequency AC techniques were 
used with currents below 10 uA. Measurements were car- 
ried out at temperatures in the range 30 mK to 3.8 K in 
the mixing chamber of a low- vibration dilution refrigera- 
tor. The superconducting magnet used was found to have 
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FIG. 3. Raw magnetometry data (change in magnetometer 
capacitance vs. magnetic field) measured at 30 mK. The lower 
left inset shows the detail at low magnetic fields. The upper 
inset is a schematic picture of the millikelvin torsion-balance 
magnetometer. The sample magnetisation causes a twist of 
the rotor which results in an imbalance in the differential ca- 
pacitor formed by the fixed and rotor plates. 



a remanent field of 26 mT, which has been subtracted 
from all the results presented. This remanent field was 
determined by analysing the (well-understood) dHvA os- 
cillations in a GaAs heterojunction, and confirmed using 
a Hall-probe measurement. 



III. RESULTS 

A. Hall Effect 

Figure H] shows the Hall resistivity p xy as a function of 
magnetic field at 30 mK. A large p xx signal, caused by 
imperfect contact geometry, was removed by subtract- 
ing the Hall resistances measured in forward and reverse 
magnetic fields. To understand the behaviour of the non- 
oscillating part of the Hall resistivity, we apply Drude 
theory, writing the conductivity as a sum of electron and 
hole contributions 
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where n e , nh and p e , ph are carrier densities and mobil- 
ities. In a non-compensated system (n e 7^ nh), the Hall 
resistivity in strong fields (p e ^B 3> 1) is determined 
by carrier density imbalance, p xy = B/e(rih — n e ). In 
graphite, however, the carrier densities are compensated, 
n e « nh, and the Hall resistivity is due to non-equal 
mobilities, p e 7^ ph- Assuming perfect compensation, 




FIG. 4. Hall resistivity p xy plotted as a function of magnetic 
field at 30 mK . Downward pointing features (marked H) occur 
when the Fermi energy passes through hole Landau levels; 
upward ones (E) occur when the Fermi energy passes through 
electron Landau levels. 



= rih = n, we obtain the Hall resistivity 
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(Interestingly, one can show that this relation holds at 
arbitrarily strong non-quantizing magnetic field.) Equa- 
tion[S]shows that the negative slope of the Hall resistivity 
in Fig. [4] is due to lower hole mobility. Using the car- 
rier density n = 3 x 10 18 cm~ 3 obtained by McClure 22 
and the average slope of p X y/B « 0.7/iOmT -1 , we find 
the relationship between electron and hole mobilities to 
be u e « 2fih- 

The oscillations in p xy may be explained as follows. 
When the Fermi energy passes through a hole (electron) 
Landau level, scattering is enhanced and /ihil^e) is re- 
duced causing a downward(upward)-pointing feature in 
p xy . Thus all the dips below 4T are due to hole Lan- 
dau levels, while the twin peaks near 7.5 T are caused 
by a spin-split electron Landau l evel. Thi s interpreta- 
tion is in agreement with that of IWoollaml and enables 
us to confirm that the carriers with the higher oscilla- 
tion frequency (higher value of Bq) are electrons while 
those with t he lower frequency are holes ( contrary to the 
assertion of lLuk'vanchuk and Kopelevichl L 



B. De Haas— van Alphen Effect 

A typical plot of magnetic susceptibility vs l/B, taken 
at 30 mK is shown in Fig. \5\ The magnetic field in this 
and subsequent figures is the component of the applied 
field perpendicular to the graphite planes. The magne- 
tometer was calibrated in situ using the known electro- 
static force between the capacitor plates^ and the raw 
capacitance data (Fig. [3]) converted to torque and hence 
to magnetic moment, using r = m x B. The magnetic 
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FIG. 5. Oscillatory magnetic susceptibility Ax as a function 
of inverse magnetic field at 30mK. The arrows indicate the 
ranges of B^ 1 used to obtain the FFTs of Fig. [6] 



moment is converted to magnetisation and then differen- 
tiated to give magnetic susceptibility. The data in the 
figure have been filtered to remove high-frequency noise. 
However, all the analysis presented in this paper has been 
carried out on unfiltered data. The fast Fourier transform 
(FFT) algorithm used later requires that the susceptibil- 
ity vs l/B data be interpolated to an integer power of 2, 
in this case 2 14 points. 



1. Magnetic-Field and Temperature Dependences 

Figure \6\ shows the FFT of the 30 mK magnetic sus- 
ceptibility data from 0.27 T" 1 to 5T" 1 . Over this field 
range, both hole oscillations (FFT peak 4.36 T, marked 
H) and electron oscillations (5.97 T, marked E) are visi- 
ble, as well as their second harmonics. The insets show 
the FFT amplitudes for the high-field ((a), from 0.27 T" 1 
to 2T- 1 ) and low-field ((b), from 2 T" 1 to 10 T" 1 ) ranges 
respectively. In the high-field FFT both peaks are of 
roughly equal height, whereas in the low-field FFT, the 
electron peak dominates. These data demonstrate that 
the (low-frequency) hole oscillation amplitude reduces 
more rapidly with decreasing field than that of the elec- 
trons. This implies either that the hole Landau levels are 
more closely spaced than the electron ones (the cyclotron 
masses satisfy the condition m/, > m e ), or that holes are 
more strongly scattered than electrons. 

In order to distinguish between these two interpreta- 
tions, we have examined the temperature dependence of 
the dHvA oscillations in Fig. [7] The oscillations at high 
l/B, which are predominantly due to electrons, diminish 
rapidly as the temperature increases, while the oscilla- 
tions at low l/B (holes) remain. This is shown explicitly 
in the FFTs of these data, in Fig. [SJ The fact that hole 
oscillations are more robust to temperature than those 
of electrons rules out the suggestion above that the hole 
Landau levels are more closely spaced than those of the 




FIG. 6. FFT of the 30 mK magnetic susceptibility (Fig. 
using data in the range from 0.27 T _1 to 5T _1 , showing the 
presence of hole oscillations with a fundamental field of 4.36 T, 
marked 'H', and electron oscillations of approximately 5.97 T, 
marked 'E'. The second harmonics of the two carriers are also 
visible at twice their fundamental frequencies. To illustrate 
the magnetic field dependences of the two sets of oscillations, 
the inset (a) shows the FFT from 0.27 T" 1 to 2T" 1 ('High 
field' on Fig. [S| where hole and electron peaks have roughly 
equal amplitude, while inset (b) shows the FFT for the data 
range from 2T" 1 to 10 T" 1 ('Low field' on Fig. O demon- 
strating the dominance of the electron oscillations at low B. 
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FIG. 7. dHvA oscillations at a range of temperatures. Oscil- 
lations at high B^ 1 are predominantly due to electrons and 
these disappear rapidly as temperature increases leaving only 
the hole oscillations at low B~ x . Traces for successive tem- 
peratures are offset for clarity. 



electrons, and indeed indicates that the opposite is the 
case, nth < rn e . The observation that the high frequency 
carriers have larger mass than the low frequency ones 
is in agreement with Soulc et al.— who found cyclotron 
masses of 0.039 mo and 0.057 mo for the low-frequency 
and high-frequency carriers, respectively (although they 
assigned the wrong carrier types to these, as discussed 
above) . Thus the only possible interpretation of Fig. [B] 
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FIG. 8. FFT plots of the dHvA oscillations using data in the 
range from 0.27 T _1 to 2T~\ for a range of temperatures. 
The high frequency (electron, E) peak is suppressed more 
rapidly with temperature than the low frequency (hole, H) 
peak. 



is that the holes experience stronger scattering. This is 
consistent with our assertion from the Hall data (Section 
MI Ap that /i e > fj, h . 

Various effects could cause different scattering rates for 
holes and electrons. The scattering rate is determined by 
the integral of the density of states over the Fermi sur- 
face and the different densities of states for electrons and 
holes, caused by the different shapes of the Fermi surface 
at the K point and near the H points, will give different 
rates. The different natures of the carriers, DF or SF, 
also affect scattering. It is known that in graphene and 
carbon nanotubes the chiral nature of the DFs suppresses 
back-scattering* 2 ^ Further experiments on different types 
of graphite would help to clarify the causes of this effect. 

We can obtain quantitative information about the car- 
rier effective masses and about scattering, from an anal- 
ysis of the damping of the oscillations in magnetisation 
with magnetic field and temperature. The amplitude of 
the magnetisation oscillations is given by: 
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where \ = ^^hsT /huj c , lo c = eB/m* and r q is the quan- 
tum lifetime. Determining r q is only possible for the 
carrier whose oscillations are dominant, in this case elec- 
trons. We find that, for electrons, m* = 0.046 ±0.003 too 
and T q = 1.7 ±0.1 x 10~ 12 s. This effective mass is some- 
what lower than the value reported by Soule et al.— If we 
assume that r q is the same as the momentum relaxation 
time (the assumption of short-range scattering) then we 
obtain a mobility /i e = 6.7 ± 0.4m 2 (Vs) _1 . 
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FIG. 10. FFT phase analysis of the 30 mK dHvA results for 
electrons - a phase of —0.25 ± 0.25 [2tt] is observed. The error 
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FIG. 9. Intercept phase analysis of the 30 mK dHvA oscil- 
lations, (a) Analysis of the oscillations of electrons. Peak 
position vs Landau-level index n is plotted up to 5T _1 . The 
intercept (see inset) indicates a phase of — 0.20±0.05 [2tt]. The 
open circles show how the low-index raw data deviate from 
the straight line, (b) Analysis of the hole oscillation extrema 
(peaks and troughs) indicating a phase of —0.48 ± 0.06 [2tt], 
consistent with 2D SFs. 



2. Carrier Natures - Phase Analysis 

From the peaks in the FFT amplitudes the funda- 
mental fields for the charge carriers have been identi- 
fied. They are: low frequency holes at 4.36 ±0.05 T 
and high frequency electrons at 5.97 ±0.01 T. Following 
the approach o f Luk'yanchuk and Kopelevich— and of 
ISchneider et al\ £ we can also extract the phase of the 
oscillations to determine the nature of the two types of 
carrier - whether they are DFs or SFs. Two types of 
analysis can be applied: a plot of the 1/B positions of 
the oscillation extrema vs Landau-level index yields an 
intercept from which the phase can be extracted; and the 
phase can be obtained directly from the FFT. 

The phase is to be obtained using Eq. ^Q which ap- 
plies to the fundamental frequency of the susceptibility 
oscillations. Since the raw data are not sinusoidal, they 



need to be filtered in order to perform an intercept anal- 
ysis. Fourier filtering was used, with a pass band centred 
approximately on the peak frequencies obtained from the 
FFT, and a band width sufficiently narrow to exclude the 
other peaks in the FFT. Various combinations of band 
width and centre were tested, as well as other filter meth- 
ods, and the intercept results were found to be fairly in- 
sensitive to these choices. 

Figure |H] shows the results of intercept phase analysis 
for the 30 mK data of Fig. [5] Low-index oscillation ex- 
trema (up to n — 4) in the raw data (open circles) are 
shown in the electron plot (a) to illustrate the deviation 
which has been commented up on by Schneider et al- 
and by ISmrcka and Goncharukl i^ but are not used in 
the fit. No such deviation is apparent in the hole data. 
An intercept of zero corresponds to a total phase of zero 
(2D DFs from Table ILT|) . while an intercept of —0.5 cor- 
responds to a total phase of — 7r (2D SFs). The results 
are summarised in Table IIIII From Fig. the intercepts 
(in units of 27r) are: for electrons, — 0.20± 0.05, not con- 
sistent with any of the assignments in Table HU for holes, 
— 0.48±0.06, consistent with their assignment as 2D SFs. 
The error bars quoted for these intercepts reflect both the 
scatter of the straight-line graphs and small changes in 
intercept dependent on filtering parameters. 

Phase information was also extracted directly from the 
FFTs. The phase was found to vary rapidly with field 
close to the FFT peaks (Fig. [TO]) , and consequently the 
uncertainty in the phase was governed by the accuracy 
with which one could determine the peak positions, which 
in turn was controlled by the range in 1/B of the suscepti- 
bility data. With the data spacing of Fig. [10] we estimate 
the uncertainty in the position of the FFT peak to be 
about ±0.1 T leading to an error in the phase of about 
±7r/2. Consequently, this gives a reasonable indication 
of carrier nature, but is not sufficiently accurate to de- 
termine dimensionality. Nevertheless the results of this 
method are in very good agreement with those of the in- 
tercept method (Table IIIII) : the phases (in units of 27r) 
are: for electrons, —0.25 ± 0.25, the error bar being too 
large to determine the nature of the carriers; for holes, 
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4.36 ±0.05 
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TABLE III. Experiment conclusions for the nature of carriers in graphite. 
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FIG. 11. (a) A(t„ as a function of inverse magnetic field at 
30 mK. The upper inset shows the raw R xx data as a function 
of magnetic field. The lower inset is the FFT of Aa xx ; the low 
frequency (hole) peak and its first harmonic are visible but 
the high frequency (electron) peak is not. (b) Intercept phase 
analysis of the SdH oscillation extrema (peaks and troughs) . 
Only holes were could be analysed and first-harmonic filtered 
data were used as discussed in the text, yielding a phase of 
—0.48 ± 0.01 [27r], consistent with 2D SF carriers. Raw ex- 
trema are also shown but exhibit systematic shifts to the right 
and to the left of the filtered data due to the presence of higher 
harmonics. 



Figure lllal shows the oscillating part of a xx , Aa xx , 
as a function of 1/B at a temperature of 30 mK. It was 
extracted from the raw longitudinal resistance R xx data 
(upper inset to Fig. Illa[) using 



Aa xx = A 



rxx * rxy 



Ap x 
rxi 



-AIL 



(7) 



The oscillations up to 4 T correspond to those seen in 
the Hall effect, which have been identified as hole os- 
cillations. Only the two features around 7T are associ- 
ated with electrons. The FFT of Aa xx (lower inset to 
Fig. Illa[) shows a clear peak due to hole oscillations but 
no evidence of electrons. This is consistent with the lack 
of electron-related features in the raw data, but contrasts 
with the analysis of the susceptibility data. 

To understand the absence of electrons in the SdH 
data, we note the following. For a conductor with only 
one carrier type, the SdH oscillations and dHvA oscilla- 
tions are both related to the oscillations in the density of 
states, and are therefore proportional to each other (see, 
e.g. Section 11.1 of Abrikosov— ) : 



6a x 



m*B, 



D 



dM 
~dB 



(8) 



-0.45 ± 0.25, consistent with their assignment as SFs. 



where m* is the cyclotron mass, and Bq is the fundamen- 
tal field. Note that the oscillating correction to the con- 
ductivity is proportional to the conductivity itself. The 
latter is inversely proportional to the mobility in strong 
fields: a xx oc l/(fj,B 2 ). Therefore, the amplitude of con- 
ductivity oscillations at high fields, when the Dingle fac- 
tor can be neglected, is higher if the mobility is lower. 
This is because, while thermodynamic dHvA oscillations 
are entirely determined by the density of states, SdH os- 
cillations are caused by the enhancement of the scattering 
rate by the peaks in the DOS. For this reason, they also 
include the scattering rate as a prefactor. 

In graphite, different carrier types contribute indepen- 
dently to conductivity oscillations. The relative mag- 
nitude of such contributions depends on the scattering 
rate for a given carrier type. Thus, high-mobility elec- 
trons should contribute less to the SdH oscillations than 
low-mobility holes. 

Figure lllbl shows the intercept phase analysis for the 
SdH data. The SdH data have significant harmonic con- 
tent in the range over which the analysis is carried out 



as can be seen from their FFT (lower inset of Fig. Illal) . 
The same filtering process as described for the dHvA data 
above was used to obtain extrema positions in Fig. Illbl 
The raw extrema are shown for comparison. A phase of 
—0.48 ±0.01 [2ir] is in excellent agreement with the anal- 
ysis of the dHvA data (Fig. Illbft once again identifying 
the holes as 2D SFs. FFT phase analysis of the SdH 
data yields a hole phase of —0.49 ± 0.25 [2tt]. This is in 
agreement with the intercept phase and the phases from 
the dHvA analysis, and indicates that holes are SFs, but 
again is subject to a substantial error. 

The amplitude Ap xx of the resistivity oscillations has 
the same magnetic field and temperature dependences as 
AM (Eq. [6J. Analysis of this amplitude gives, for holes, 
to* = 0.033 ± 0.003 too and r q = 2.1 ± 0.2 x 10~ 13 s. 
This effective mass is in reasonable agreement with the 
value reported by Soule et al. — Assuming that t ? is the 
same as the momentum relaxation time we obtain a hole 
mobility fi h = 1.2 ± 0.1m 2 /Vs. The dHvA amplitude 
calculation for electrons f Section Till B ip gave a mobility 
of fi e = 6.7 ± 0.4m 2 /Vs. The amplitude analyses there- 
fore confirm qualitatively the conclusion from the Hall 
effect data (Section IIII A[) that /i e > fih- Quantitative 
agreement is not expected because the Hall effect is sen- 
sitive to relaxation times averaged over the entire Fermi 
surface while magneto-oscillatory effects primarily probe 
behaviour at the extremal cross-sections. 



IV. CONCLUSIONS 

We have carried out Hall, SdH and dHvA measure- 
ments on the same sample of HOPG, at temperatures 
down to 30 mK. Our Hall effect experiments confirm the 
carrier types established by Schroeder et al.—: holes near 
the H point of the graphite Brillouin zone are responsi- 



ble for the low- frequency component of the dHvA oscilla- 
tions, while electrons at K give rise to the high-frequency 
component. Analyses of the Hall, dHvA and SdH effects 
all indicate that the holes are subject to stronger scatter- 
ing than the electrons. The origin of this observation will 
be the subject of further investigations of various types 
and grades of graphite. The temperature dependence 
of the dHvA effect indicates that the holes have larger 
Landau-level separations than the electrons, a conclu- 
sion confirmed by amplitude analyses of dHvA and SdH 
oscillations. 

The analysis of the phase of the dHvA and SdH os- 
cillations suggests that the low-frequency carriers (holes) 
are SFs, in agreement with the previous studies sum- 
marised in Table U and Fig. [2] The phases are consis- 
tent with the holes being two-dimensional, in agreement 
with the HOPG experiments of Schneider et al.— but 
not with other experiments. This probably reflects the 
fact that the interaction between atomic layers is differ- 
ent in the different samples used. The situation with 
the high-frequency carriers (electrons) is more compli- 
cated. Our results show an intermediate value of the 
phase, consistent with neither DF nor SF. Other au- 
thors ( Williamson et ai\ ^ Schneider et al.—) have also 
reported an intermediate phases for the high-frequency 
carrier. It is interesting to note that Mikitik and Shar- 
lai 5 pointed out that intermediate phases are possible 
for minority holes in situations close to magnetic break- 
down, though this could not occur for majority carriers 
with standard values of the SW parameters. Intermedi- 
ate phases might also be the result of sample inhomo- 
geneity, in which electrons are SF and DF in different 
parts of the same sample. It is not clear why this should 
affect electrons more than holes though. Clearly graphite 
remains an interesting and intriguing material and more 
measurements on a range of graphite samples are needed 
in order to determine the nature of its electrons. 
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